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ABSTRACT 


The propagation of fire generated smoke into a shipboard space with a geometric 
interference has been modeled using commercial software from the Computational Fluid 
Dynamics Research Corporation (CFDRC). This study was based on the dimensions of 
compartment 01-163-2-L and the installed ladder aboard an Arleigh Burke Class Flight 
IIIA Destroyer. A test model was run which validated the hindrance of fluid flow by a 
geometric interference. Smoke propagation scenarios were run in the shipboard 
compartment model. The results of the first scenario showed that smoke propagation is 
limited by the geometric interference. The results of the second scenario showed that 
smoke that 1s directed vertically 1s diverted by the geometric interference. The overall 
goal of this study is to show that computational fluid dynamics software can successfully 


model smoke propagation in shipboard spaces with a geometric interference. 
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I. INTRODUCTION 


A. BACKGROUND 


Throughout US naval history, fighting a shipboard fire has been the number one 
damage control priority. Aboard a ship, evacuation from a fire is not a choice. In a 
recent event, the USS STARK (FFG 31) was struck by two Iraqi Exocet missiles while 
on patrol in the Arabian Gulf. Both missiles entered on the port side of the STARK but 
most of the damage was on the starboard side [Figure 1]. The first missile failed to 
detonate but spread deadly burning propellant in its path. The burning propellant 
generated extremely high temperatures causing thermal damage and enormous amounts 
of smoke. The second missile detonated within the skin of the ship leaving a gaping hole 
in the hull. This hole fed oxygen to the fires caused by the extreme temperatures. Decks 
and electrical cableways melted from the 3000° F(1922° K) temperature produced from 
the burning propellant. The crew of the STARK fought the high temperature fire that 
produced tremendous amount of smoke and toxic fumes. 

‘The heat and smoke were tremendous,’ LT Carl S. Barbour 
recalled. For example, when I cracked the hatch from the mess decks by 

the scullery, it felt as if the fire was right there. Yet, we didn't find flames 

until we got all the way to the rear of the CPO berthing.[Ref 1] 

The crew’s fire fighting capabilities were limited by the dense smoke, toxic fumes 
and the US Navy’s fire fighting technology. It was impossible to fight the fire through 
the smoke and fumes as the missile propellant burned unabated. Fighting this weapon- 


induced fire was new to the US Navy. Today, weapon-induced fires remains a hazard 


even with the newest fire fighting technology aboard ship. 
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Figure 1 from Ref 2. USS STARK (FFG 31). Starboard Side Damage. 





B. DESCRIPTION OF THE PROBLEM 


Modern technology today, allows computer simulation to enhance engineering 
research and development. In U.S. Naval applications, computer simulation modeling of 
fire and smoke propagation in current shipboard spaces would facilitate plans of action in 
damage control, and future ships could be designed space by space for maximum safety 
in fire and smoke protection. These computer modeled ship spaces could designate 
where fire fighting and life saving equipment should be placed for easy access in case of 
an emergency. 

Most spaces aboard a ship are designed for maximum use of the equipment and 
machinery that are designated to that space. Equipment and machinery in a space 
become obstacles or fuel for a fire during a blaze where heavy and toxic smoke is 
present. During a damage control assessment of how to attack a fire, knowing how the 
smoke will propagate in a space will enhance the Damage Control Assistant's decisions 
of how his repair parties will combat a fire. As desktop personal computer speed and 
memory increases every year, the capability of predicting smoke propagation for every 
space in a ship is not far off in the future. 

Each new class of ship becomes more technologically advanced, thus less 
personnel are needed to man it. It is planned that only 95 personnel will safely and 
efficiently operate the Surface Combatant of the 21“ Century (DD 21). 

Improvements in design can contribute significantly in obvious 

ways by reducing both the susceptibility and the survivability of platforms 

before an attack occurs. However, the benefits of applying technology to 

the reactive effort following a successful enemy attack, are hard to 

quantify. ... there are a number of damage control experts who doubt that 


technology can contribute to DC so comprehensively. They argue the 
impossibility of being able to predict the location or degree of damage and 
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that the path to increased survivability lies in the direction of adequate 
manpower and better platform design. [Ref 3] 


With a small number of the crew to man fire parties, it is essential that 
knowledge of fire, heat, and smoke propagation in ship spaces is readily available to all 
personnel. The results of this study could offer future design engineers the data on how 
fire generated smoke propagates in shipboard space geometry and how the solid 
interferences within the space will affect it. 


C. PREVIOUS WORK 


Jones and Walton [Ref 4] took their knowledge from their study of fire and smoke 
propagation in buildings and applied them to ships. Their methodologies were developed 
for civilian structures, so algorithms for stairways and elevators were changed for 
hatches, scuttles, and ladders. At the David W. Taylor Naval Ship Research and 
Development Center (DTNSRDC), they simulated a scenario where a 1-megawatt fire 
was located on the front starboard locker and berthing space. Zone modeling was used for 
this study. The simulated fire was caused by a missile hit and the burning of unspent 
solid propellant. These simulations were run in 1983 and were evaluated with limited 
processor speed and software. They determined that model simulations that had been 
developed for predicting fire and smoke propagation in buildings were similar in ships. 

In 1985 Jones [Ref 5] studied fire and smoke propagation in multi-compartmental 
spaces implementing 2 computer programs. The BUILD software program was used to 
generate the model configuration, and the FAST software program was used to run fire 
scenarios in the model. Scenarios were run to emulate previous experiments of actual 


fires in multi-compartmental spaces. The results from the FAST program were to be 


compared to the experimental data. In previous research, Jones noted that two or more 
layers of gases formed in a compartment. He again used zone modeling and modeled 
each compartment with composition of gas layer control volumes. 

In this context, the predictive equations for the gas layers in each 
compartment result in from conservation of mass, momentum, and energy 
together with an equation of state for each compartment.[Ref 5] 

Again limited by processor speed and memory size, Jones found disparity 
between actual experimental data and computer simulated data. After running 
simulations, he found his predicted temperatures were too high and gas layer depths were 
too small. 

In 1992, Forney and Jones [Ref 6] with faster processor speed and improved 
software program CFAST, improved on modeling smoke movement through 
compartmented spaces. They successfully presented the radiative and convective heat 
balance terms which affected smoke flow through buildings. Their work emphasized the 
movement in a space of toxic gases that are generated in a fire. Their predictions of 
radiative and convective heat balances were favorable with experimental data. 

In 1993, Forney and Jones [Ref 7] further improved their smoke transport model 
from previous work. Using the CFAST software program and a faster processor, they 
modeled the movement of toxic gases from the space of origin to a distant compartment. 
They also studied smoke transport with vertical flow and with mechanical ventilation. 
Refining the radiation transport scheme which affected energy distribution and buoyancy 
forces, their improved model generated data consistent with their experimental data. 

Tatem and Williams [Ref 8] used the software program FAST and modeled 


missile propellant fires in shipboard compartments. They conducted a series of 


experimental tests of burning propellant in a steel mock-up of shipboard compartments at 
China Lake. An algorithm for the burning rate of the missile propellant was developed in 
FAST, and after each experimental test, they ran their computer simulated model. The 
China Lake simulated test results underpredicted peak temperatures and overpredicted 
heat fluxes. In a second experimental test series, missile propellant was ignited aboard an 
ex-LEANDER Class Royal Navy frigate. In these series of simulated test runs, the 
predicted peak temperatures were in agreement with the experimental data, but again heat 
fluxes were Overpredicted. 

Mehls [Ref 9] used a commercial code CFD-ACE generated by Computational 
Fluid Dynamics Research Corporation (CFDRC). He modeled smoke propagation in a 
compartment aboard an Arleigh Burke (DDG 51) Class destroyer. Using his model and 
scenarios run, he was able to predict the temperatures of the mixture of smoke and cool 
air and how smoke propagates within a shipboard compartment. His model did not 


include any geometric interferences. 


D. OBJECTIVES 


The purpose of this study is to develop and examine a computationally generated 
model that can predict how smoke travels in shipboard spaces that contain geometric 
interferences. The model will be generated in the form of compartment 01-163-2-L, 
aboard an Arleigh Burke (DDG-51) class destroyer. The geometrical interference will be 
in form of a ladder with a ladderback installed. Physical properties associated with the 


smoke such as density will be simulated. 


Il. COMPUTATIONAL FLUID DYNAMICS 


A. OVERVIEW 

Computational Fluid Dynamics (CFD) software is an invaluable tool for design 
Optimization and rapid virtual prototyping for fluid transport problems. CFD computer 
simulations eliminate “trial and error" engineering and hasten the development of the 
fluid transport design and application. The commercial software package that was used 
in this research is called Computational Fluid Dynamics-Advanced Computing 
Environment (CFD-ACE) Version 6.2 [Ref 10-13] and was developed by the 
Computational Fluid Dynamics Research Corporation (CFDRC). 

CFD-ACE is an integrated package comprised of three separate, yet interactive 
codes to solve the fluid transport problem [Figure 2]. The three codes are GEOM, GUI 
and VIEW. CFD-GEOM is the processor where a model can be created from scratch, or 
the model can be imported from another CAD program. CFD-GEOM offers 
comprehensive mesh generation, enabling the user to generate structured, unstructured, 
and mixed element meshes to represent the structure of the fluid transport problem. 
CFD-GUI (Graphic User Interface) is the solver for the package. The CFD-GEOM 
meshed model is imported into CFD-GUI. In CFD-GUI, scenarios for the fluid transport 
problem are created. Scenario parameters are set by the user. After fluid properties, 
initial and boundary conditions, and interaction of species (heat transfer, turbulence, 
mixing) are set, the user designates the number of iterations to run. The conservation 
equation solutions for the model are affected by the chosen differencing scheme and by 
varying the amount of relaxation and constraints. During the run, CFD-GUI solves the 


series of equations for all the inputted parameters. After the run in CFD-GUI, the solved 
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data is imported into CFD-VIEW. CFD-VIEW graphically illustrated the results in 2-D 
or 3-D. 


B. FINITE VOLUME ANALYSIS 


The first step in the CFD analysis is to construct a geometric model over which the 
relevant fluid transport equations can be numerically integrated. The model creation is 
called domain modeling. CFD-ACE employs a structured, multi-domain, body fitted 
coordinated system approach which enables the user to simulate flows in complex 
geometric configurations. The fluid flow domain is gridded and divided into a number of 
cells known as control volumes. 

A control volume is similar to a cube with six faces and six direct neighbors 
[Figure 3]. In the finite volume approach, discretized equations are formulated by 
evaluating and integrating fluxes across the faces of each control volume. This satisfies 
the relevant conservation equations. Dependent variables are solved for at the center of 
the control volume. The values obtained are considered to prevail over the entire control 
volume. Differencing schemes of varying accuracy can be used when evaluating 
convective fluxes over the control volume. These schemes can be independently selected 
for each fluid transport variable to be solved. 

In CFD-ACE, fluid flows are simulated by numerically solving partial differential 
equations (PDE's) that govern the fluid transport variables. The mass, momentum, 
energy, turbulence quantities, mixture fractions, species concentrations, and radiative 
heat fluxes that will be solved will depend on the nature of the flow problem. The PDE’s 


are discretized on a computational grid. A set of algebraic equations are formed and 


solved. This numerical method yields a discrete solution of the flow field. The flow field 
is comprised of the values of the fluid transport variables at the grid points. 

CFD-ACE uses an iterative solution method where equation sets for each fluid 
transport variable are solved in sequence until a converged solution is obtained. In CFD- 
ACE, the SIMPLEC algorithm is used [Figure 4]. The user implements the number of 
iterations (NITER) and in the case of transient simulation, the number of continuity 
iterations (C_ ITER) to be run. NITER and C_ITER are dictated by the overall residual 
reduction obtained. At each iteration, CFD-ACE will calculate a residual for each fluid 
transport variable for all control volume cells. A reduction of five orders of magnitude in 


the residuals is needed before convergence is accepted. 
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Figure 3 from [Ref. 10]. Figure 4 from [Ref. 10]. 
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Iii. MODEL 


A. GEOMETRY 


The creation and simulations of the model were carried out using a Micron Client 
Pro Desktop computer, with 384 megabytes of RAM and an internal 12 gigabyte hard 
drive. The software used was CFD-ACE+ version 6.2, which was last updated in July 
2000. 

A test box model of a shipboard ladder with an installed ladderback was modeled 
inside an 8m (length) by 8m (depth) by 2.29m (height) box. The dimensions of the 
ladder were identical to shipboard specifications. Figure 5 is a skeletal view of the ladder 
inside the test box. Due to the simple locations of the ladder and watertight doors in the 
test box, a structured grid was used on the model. Using the structured grid, 14 control 
volumes were created. The ladder was comprised of three solid control volumes, and the 
remainder 11 control volumes were made of air. 

A model of passageway 01-163-2-L aboard an Arleigh Burke Class Destroyer 
(DDG-21) was designed. The dimensions of the model and ladder identically match the 
actual compartment. A shipboard plan view of the compartment from a ship's drawing is 
shown in Figure 6. This compartment has a variety of openings for smoke intrusion and 
will allow an assortment of smoke propagation scenarios to be studied. CFD-ACE has 
the capability for any entrance or exit in the model, when not part of the scenario, to be 
designated as a wall. Therefore, the non-activated entrances and exits have no effect on 
the results of the scenario. Using the exact dimensions and locations of the ladder, 
hatches, and watertight doors, a problem arose with the creation and orientation of the 


control volumes in the space. To remedy the situation, an unstructured grid was used on 


Mm 


the compartment. The unstructured grid allowed the ladder to be made up three solid 


contro] volumes and the entire space to be made up of a control volume of air. 
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Figure 5. Plan View of Space 01-163-2-L | 


B. GRID GENERATION 





In CFD-GEOM, a succession of steps must be followed when generating a | 
structured or an unstructured grid on a model. For both types of grids a skeletal model | 
must be created as shown in Figure 7-8. Both rough skeletal models are made up of | 


geometric lines in the line generation tools in CFD-GEOM. Both skeletal models show 





the large rectangular outer structure made up of the space, ladder, hatches, and watertight 
doors. 

For structured grid generation, edges must be placed on the lines that make up the 
model. CFD-GEOM will place a number of grid points inputted by the user on these 
created edges. A structured grid face is then created from a closed polygon that is made 
up from the selection of four edges. Each edge opposite of another edge must have the 
same number of grid points in order for the face to be made. For a circle, edges are 
represented by 4 arcs with an equal number of grid points for each arc. These arcs can 
then be picked to make up a grid face. 

For an unstructured grid generation, edges do not have to be created and put on 
the lines that make up the model. A rough surface is first placed on the polygon face 
where an unstructured grid is needed. A trimming loop is then placed on the outline of 
the polygon face. The user then must select the created rough surface to be trimmed and 
then select the trimming loop. CFD-GEOM then will trim the rough surface in the shape 
of the loop that was placed on the polygon face. This is a trimmed loop surface. The user 
then must create trimmed loop surfaces on each polygon face of the model until the entire 
model is made up of trimmed loop surfaces. After the model is made up of trimmed loop 
surfaces, a closed surface set can be made for the model. The user selects the ‘create 
closed surface set' icon in the topology section in CFD-GEOM. — Each trimmed loop 
surface must be selected. After all trimmed loop surfaces are selected, the user must 
input the selection. The closed surface set is then created. To create the unstructured grid 
surface, the ‘triangular grid' icon must be selected from the grid section in CFD-GEOM. 
The user then will select and input the created closed surface set. An unstructured grid 
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surface is then generated on the model. Rough surfaces and trimming loops can easily 


put on circles to make trimmed loop surfaces to be picked to generate unstructured grids. 
C. VOLUME GENERATION 


In a structured grid, 6 structured grid faces must be created to make up a volume 
block. To make a volume block, the user must select the ‘create a block' icon in the grid 
section in CFD-GEOM. The user will be prompted to select 6 grid faces. Once the faces 
are selected and inputted, the volume block is created. Various volume boxes make up a 
volume of a model. Opposite faces must be have the same number of grid points in order 
for the structured volume block to be generated. 

In an unstructured grid, volume cells for the model are created by first selecting 
the ‘create volume set’ icon in the topology section in CFD-GEOM. The user will be 
prompted to pick the closed surface set that were made while creating the unstructured 
grid. Once the closed surface set is selected and inputted, the volume set for the model is 
generated. The user then must then create tetrahedral cells for the volume. The 
tetrahedral icon in the grid section 1s then selected. The user will be prompted to select 
the volume set. Once the volume set has been selected and inputted, the volume cells in 
the model are created. 

D. MODEL CONFIGURATION 

The test box model shown in Figure 9 was made up of 98,224 structured grid 

cells, 50 grid faces, and 14 volume boxes. The compartment model shown in Figure 10 


was made up of 298,424 unstructured tetrahedral volume cells. 
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Figure 8. Test model with structured grids. 
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Figure 9. Model of compartment 01 
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IV. RESULTS 


The objective of this research was to use the CFD-ACE program to see how a 
geometric interference in a shipboard compartment affects smoke propagation. Three 
Scenarios were run in the generated model. Each scenario was run in steady state. All 
inputs for each scenario are shown in appendices. 

Scenario A was used to qualitatively evaluate how geometric interferences 
designed in a space can modify fluid flow. 500K air was set to enter the space through 
the forward door and exit the aft door. Figures 10 and11 display isotherms being diverted 
about the ladder. Figure 12 shows the diverted isotherm engulfing the ladder. The ladder 
has acted as a barrier and diverted the flow. 

With the success of scenario A, scenario B was run to compare the results of 
Mehls’[Ref. 9] scenario A. Figures 13 and 14 display the isosurface smoke 
concentrations of 88% and 77% respectively. In Figure 14, the isosurface begins to be 
diverted up the ladder. Figure 15 and 16 compares isosurface smoke concentrations of 
54% for scenario B and Mehls’ scenario A. The scenario B isosurface has not propagated 
as far in the space as the isosurface in Mehls’ Scenario B due to the ladder diverting the 
smoke. Figure 17, displays the isosurface smoke concentration of 40%. This isosurface 
has engulfed the ladder and is being diverted upward. The designed geometric 
interference has limited the propagation of smoke. 

Scenario C was used to study the effects the ladder has on smoke that is entering 
the space vertically. The smoke enters the hatch located on the deck and exits the hatch 
located on the overhead. Figure 18 displays smoke entering the space and the isosurface 


is 99% smoke concentration. Figure 19 displays the isosurface smoke concentration of 
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80%. The back of the ladder has diverted the isosurface around it. In Figure 20, the 
isosurface smoke concentration of 54% engulfs the ladder and is again diverted upward. 
Smoke that enters vertically is also impeded by the designed geometric interference. 

For each scenario, residual outputs decreased a magnitude of five orders. 


According to CFDRC’s criteria of convergence, scenario results were validated. 








V. CONCLUSIONS 


This study's results successfully modeled smoke propagation in a shipboard 
compartment with geometric interferences. Scenario accuracy results were validated and 
verified by the residual outputs. 

This study verified that smoke propagation in a compartment is affected by the 
geometric interference. The ladder diverted the smoke thus slowing smoke propagation 


within the compartment. 
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VI. RECOMMENDATIONS 


The following recommendations are made in continuation of this study: 

- Model a heat source in the compartment and use CFD-ACE to analyze fire 
and smoke scenarios. 

- Add adjoining compartment to analyze how fire and smoke propagation are 
affected by the ladder. 

- Mode] the compartment with more complex geometric interferences in a CAD 
software program (e.g. IDEAS) and then import it into GEOM. This will 
allow for more complex geometries to be gridded with unstructured grids. 

- Use transient time step calculations to calculate wall temperatures and rate of 


smoke propagation. 
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APPENDIX A. 


Scenario A was developed as a test scenario to see if the designed ladder with a 
ladderback would effect fluid flow. The ladder was designed inside a 8m (length) by 8m 
(depth) by 2.29m (height) box. The front watertight door is the inlet and the back 


watertight door is the outlet. Refer to the next page for the required inputs. 
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Table 1. Input data for Scenario A. 
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Figure 10. Isotherm first engaging ladder inTest Box. 
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APPENDIX B. 


Scenario B was developed as a scenario comparable to Mehls [Ref 9] scenario A. 
The inlet is located at the front watertight door. The upper half of the door is designated 
as the smoke inlet. The lower half of the door is exclusively air. Refer to the next page 


for the required inputs. 
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Table 2. Input Data for Scenario B. 
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Table 2 Cont. Input Data for Scenario B. 
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Figure 13. 88% Smoke Concentration in compartment 01-163-2-L. 
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Figure 14. 72% Smoke Concentration in compartment 01-163-2-L. 
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Figure 16. From Mehls [Ref.9] 54% Smoke Concentration in compartment 01-163-2-L. 
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Figure 17. 40% Smoke Concentration in compartment 01-163-2-L. 
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APPENDIX C. 
Scenario C was developed to study how the propagation of smoke in the vertical 
direction is affected by a geometric interference. The bottom scuttle is designated a 
smoke inlet while the top scuttle is designated a smoke outlet. All the watertight doors 


are designated as walls. Refer to the next page for the required inputs. 
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Table 3. Input Data for Scenario C. 
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Table 3 Cont. Input Data for Scenario C. 
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Figure 18. 99% Smoke concentration in compartment 01-163-2-L. 
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Figure 19. 80% Smoke concentration in compartment 01-163-2-L. 
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Figure 20. 54% Smoke concentration in compartment 01-163-2-L. 
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